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Abstract 

We study the antiferromagnetic O(A^) model in the F4 lattice. 
Monte Carlo simulations are applied for investigating the behavior of 
the transition for = 2, 3. The numerical results show a first order 
nature but with a large correlation length. The — > 00 limit is also 
considered with analytical methods. 
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1 Introduction 



The antiferromagnetic formulations of field theories in four dimensions have 
been recently paid considerable attention [||, H, ||, ||. The hope is to give 
some insight into the well known triviality problem in field theory 0. Also 
there are other interesting phenomena as the apparition of new particles § . 

A spin model, in a simple cubic lattice with first neighbor interactions, 
becomes antiferromagnetic if the coupling is negative. However, with some 
exceptions 0, , a simple staggered transformation maps the antiferromag- 
netic phase into the usual ferromagnetic one. 

To obtain a non-equivalent antiferromagnetic phase one has to include 
further couplings or modify the lattice geometry (see for instance ref 

Perhaps the simplest method to obtain non-trivial antiferromagnetism in 
four dimensions is to work in an F4 lattice. It is defined by taking out the 
odd sites (the sum of the coordinates is odd) of a simple hypercubic lattice. 

Four dimensional antiferromagnetic O(A^) models have been already stud- 
ied by Monte Carlo(MC) means in this lattice. The 0(1) model (Ising model) 
was considered in reference 0: a weak first order transition was found. A 
study of the 0(4) model appears in reference 0: in the range of the lattice 
sizes simulated, the behavior pointed to a second order transition. 

In this letter we consider the intermediate cases: 0(2) and 0(3), to know 
if the order of the transition changes with N . We will give evidence that the 
transitions are in both cases first order, but the numerical difficulties grow 
with N . We also present an analytical study in the N ^ 00 limit. 



We label the coordinates of an F4 lattice as a set of integers {x, y, z, t} such 
that x + y + z + t is even. We consider the action 



where the sum runs over 24 pairs of nearest neighbors and the field, is 
a normalized component real vector. We work in a hypercubic lattice of 
size V = L^/2 with periodic boundary conditions. 

In the ferromagnetic region {(3 > 0) this model is expected to belong to 
the same universality class of the simple cubic model: it presents a second 



2 The Model 




(1) 



<i,j> 
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order transition with mean field critical exponents. 

In the antiferromagnetic sector {(3 < 0) the system also presents an or- 
dering phase transition, but the structure of the ordered phase is much more 
complex. In fact, it can be easily checked that the ground state presents 
frustration. Moreover, the ordered vacuum is not isotropic. 

The independent order parameters we can construct with periodicity 2 

are 

x,y,z,t 



x,y,z,t 



(2) 



V 

x,y,z,t 



where the sums are extended to all the F4 sites. The dots stand for the other 
3 combinations of hyperplanes and 2 of planes. 

We label the site / = 1, . . . , 8 inside the 2"^ elementary cell with its Carte- 
sian coordinates X,Y,Z,T = 0,1. In practice we measure the 8 different 
magnetizations associated to a given position in the elementary cells, rrii, 
defined as the normalized sum of the magnetization for all 2^ cells for each 
of the 8 sites. 

x,y,z,t 
x — X, . . . even 

The quantities (§) can be expressed as linear combinations of these mag- 
netizations. 

The mean magnetization and the susceptibility are defined respectively 

as: 

The Binder Cumulants are defined in a such way that Vm|^=o ~ ^ 
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For the connected susceptibility we use the definition 



1=1 



3 Critical behavior 

For an operator O that diverges as {(3 — /3c)~^°, its mean value at a coupling 
/3 in a size L lattice can be written, in the critical region, assuming the 
finite-size scaling ansatz as 0] 

0{L, P) = L^ol^ {Fom, P)/L) + 0(L-)) , (7) 

where Fq is a smooth scaling function and uj is the universal leading corrections- 
to-scaling exponent. In order to eliminate the unknown Fq function we can 
measure at the coupling where Fq presents a maximum as for Xc or for the 
specific heat Cy- Another method is to study the behavior of the quan- 
tities 

Qo = 0{sL,P)/0{L,(3). (8) 
We use VmL as scaling variable |^. It is direct to obtain 

Qo\q.,,.=s = s''-'' + 0{L--). (9) 



We use that = 7 and xa^iog(M) = 1 to obtain the critical exponents. 

The expected FSS behavior of a first-order transition |10] corresponds to 
apparent exponents: u = a = 1, 7=1. 



4 The Simulation 

We will consider in this letter the cases = 2, 3. We have used a Metropolis 
algorithm followed by Nq overrelaxation steps as update method. No depend- 
ing on the model and lattice size. We have worked in lattice sizes up to 48. 
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L 


/ 


No 




#ofr 


/5 


4 


20 


3 


0.775(9) 


100000 


-0.352 


6 


20 


3 


1.144(14) 


87000 


-0.352 


8 


20 


3 


1.63(4) 


55000 


-0.352 


12 


20 


3 


3.16(6) 


60000 


-0.352 


16 


20 


3 


6.5(5) 


11400 


-0.352 


24 


25 


4 


15.8(7) 


5000 


-0.3516 


32 


24 


7 


32(1) 


1700 


-0.3513 


48 


32 


7 


85(7) 


312 


-0.35125 


4 


20 


3 


0.871(11) 


98000 


-0.53 


6 


20 


3 


1.19(2) 


82000 


-0.53 


8 


20 


3 


1.52(3) 


98700 


-0.53 


12 


20 


3 


2.07(4) 


63000 


-0.53 


16 


20 


3 


3.15(6) 


48000 


-0.53 


24 


24 


5 


4.99(14) 


29000 


-0.5287 


32 


28 


6 


8.7(4) 


5000 


-0.5287 


48 


28 


6 


21(1) 


2000 


-0.5286 



Table 1: Description of the simulation in the critical region for 0(2) (upper 
part) and 0(3) (lower part). We report the lattice size, the frequency of 
measures, number of over-relaxation steps for each Metropolis one, autocor- 
relation time (r) for iterations in r units, and the coupling. 



We used for the computation the dedicated machine RTNN, consisting in 
16 Dual Pentium Pro units. For the largest lattices we parallelized, using 
shared memory, in each dual motherboard. Every / sweeps we measured the 
energy, the specific heat and the 8 period-two magnetizations m/. In table || 
we report the parameters of the simulation at the critical region. 

In order to extrapolate to the neighborhood of the critical point, we used 



the usual reweighting method |TT[. We also simulated at the maxima of the 
specific heat because the region where the extrapolation was reliable was not 
large enough. These simulations were about one fourth of the total CPU 
time. 

The errors were computed with a jack-knife method, performing 50 blocks 
in order to achieve statistical error bars within 10% of precision. 
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5 The Vacuum 



For large |/?|, all the magnetizations mj go to unitary vectors, confirming the 
assumption of a period 2 vacuum, so that we can restrict our analysis only 
to a unit cell. We also find that A^afp is non zero, while Mp and M^fr 
vanish. 

When Mf = Mafk = , it follows that mxYZT = 'n^i-x,\-Y,i-z,i-T- 
We checked that the cosine between tuxyzt and mi_x,i-y,i-z,i-T goes to 
1, with the expected behavior in the broken phase corresponds to AFP 
order. 

This ordering is also found near the transition. So we can restrict our- 
selves to study only 4 independent spins. 

At T = the frustrated ground state is described in it consists in two 
couples of antiparallel spins, but the angle between couples is free. The differ- 
ent choices for the couples determine the plane for AFP symmetry breaking. 

At T > we do not know whether there is a privileged angle or not. It 
is necessary to determine the pattern of the symmetry breaking. 

A very important point is to determine if the two couples are aligned or 
not. To clarify this point we construct the following tensor: 



where the superindices run for the components of the N vectors and the sum 
over the four independent magnetizations in the elementary cell. It is clear 
that if the two couples are aligned, the four tensors as well as the sum tensor 
can be simultaneously diagonalized, and so we expect in the broken phase a 
non-zero value for the largest eigenvalue and zero values (up to effects) 
for the rest {N — 1) of them. 

If this holds true, we also expect in the critical region a L^'^^^'^ behavior 
for the biggest one and a for the others. This will be checked in the next 
section. 

In table we show the eigenvalues for both models in the broken phase. 
We note that the largest one goes to a non-zero value while the rest go to 
zero. The fit parameters are obtained with a linear extrapolation in for 
the former case and in L""^ for the latter. 




(10) 



K=l 
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Eigenvalue 


xVd.o.f. 


Value 


-^max 


0.16 


0.20346(15) 


^min 


0.13 


-0.0(2.2) xlO-'^ 


-^max 


0.11 


0.18986(7) 


■^med 


0.89 


2.3(2.7)xl0-^ 




0.05 


-1.7(1.2)xl0-^ 



Table 2: Eigenvalues for the tensor ([T0|) for the 0(2) model (upper part) at 
/3 = —0.37 and for the 0(3) model (lower part) at /3 = —0.57. 

6 Critical exponents 

A determination of the critical exponents is obtained by studying the height 
of the peaks of Cy and Xc- For a first order transition both quantities should 
diverge as the volume {a/u = 4, 7/z/ = 4). For small lattices it is usual to 
find the apparent critical exponents of a weak first order transition: ajv = 
1,7/z/ = 1 (see ref. [0), which are precursors of a first order transition. 

To analyze the divergence of Cy a bilogarithmic plot is not adequate due 
to the presence of a non- negligible constant term. In order to compare with 
the first order behavior it is better to plot Cy and function of several 

powers of L. This is done in fig. 1. We remark that in the 0(2) case for L in 
the interval [8,24] ( [12,32] for 0(3)) there is an excellent linear fit for n = 1 
which is the value predicted in a weak first order transition. This frequently 
produces a misunderstanding of the order of the transition. However, it is 
clear from fig. [l| that this is a transient effect: for the larger lattices the 
divergences are faster than linear, and presumably they would reach the first 
order behavior for very large lattices. 

The susceptibility also shows a fast divergence. Although we are not able 
to observe the asymptotic first order behavior the trend seems rather clear. 

A more accurate measure of the critical exponents can be obtained from 
eq. (|^). We have always used the ratio s = 2. In fig. |^ we plot several 
determinations of exponents using different operators. We remark that there 
is a systematic error in the ajv determination because of the analytic term 
in Cy. We observe no asymptotic behavior in all cases although the values 
in the larger lattices are hardly compatible with a second order transition. 

In the upper right part of fig. ^ we plot the exponents related with each of 
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0(2) 0(3) 

- 1 I I I I I I I I I I I I I I I I I I I I I I I I I u 2.0 I — \ — \ \ \ \ \ \ \ \ \ \ \ \ \ I 




(-^'/-^max)" (-^-/-^max) 



Figure 1: Specific heat and connected susceptibility for 0(2) and 0(3) as a 
function of several powers of the lattice sizes. 



the eigenvalues of the matrix (|iy). While the maximum eigenvalues should 
behave as (with /?/z/ = for a first order transition), the others should 

go to zero as L~'^. The latter can be used as a control of when the asymptotic 
regime is reached. We observe that we are far from this regime but the first 
order limit seems rather clear. 

A comparison between the curves for 0(2) and 0(3) shows a roughly 
similar shape, differing in a horizontal shift that corresponds to multiplying 
the lattice size by a factor near 2. This fact can be understood as a correlation 
length at the critical point which is twice larger for 0(3) than for 0(2). 
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Figure 2: Critical exponents for 0(2) (solid lines) and 0(3) (dashed ones) 
measured through the relation (^). The filled circles (diamond for the non- 
maximum eigenvalues) mark the first order limit. 

7 Critical point and energy histograms 

To calculate an estimation for the critical coupling we study the (3l values 
where Vm{2L, Pl) = Vm{L,Pl)- In both models, these Pl for the largest 
lattices are compatible within the error bars. We can fit to the functional 
form — /5c (oo) oc in order to estimate the error bars. We perform 
these fits with the full covariance matrix. In both cases we obtain fitting for 
L > 6 a wide valid range for x and very good (x^/d-o.f. = 1.9/2 for 0(2) 
and x^/d.o.f. = 0.8/2 for 0(3)). The results are compatible with the values 
for the largest lattices. We get 

/3^(oo)°(^^ = -0.351216(10) (11) 
/3^(oo)°(^^ = -0.52857(2). (12) 
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Figure 3: Normalized energy histogram at /?c(oo) for both models. 



Finally let us comment on the energy distribution of the configurations. 
A direct check of the first order character of a transition is the observation of 
a latent heat. Unfortunately, a sharp double peak structure can be observed 
only when the lattice size is much larger than the correlation length at the 
critical point. In figure ^ we show the energy histograms for both models 
at Pc- In the 0(2) case we note that the width of the energy distribution is 
nearly constant for the larger lattices, being an indication of the existence 
of a two peak distribution that cannot be resolved. In the 0(3) case, up to 
L = 48 there is not a similar behavior. 
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8 The N ^ oo limit 

The partition function of the model can be written as 

Z = zJ _ 1) ^-NP'H 

^ j 

= /'n^^^.^e^^^^^'"^'''^^^'^*^'^"''''^ (13) 
J . 2n 

where /?' = P/N, and Zq is a normahzation factor such that Z ^ 1 when 
P ^ 0. We have introduced the conjugate parameters (yj,Xj to give an 
integral representation of the constraint = 1 [[1^ . 

Writing the quadratic form in the exponent of eq. ( [T^ ) as 

-Y^^nQnm(pm, (14) 

n,m 

the integration over $ yields 

Z = Zol^] / n^«.eME.2(A,+ia,)-THnQ}_ ^^5) 

In the limit oo, a variational equation with respect to 2{Xk + iak) 

gives 

1 = (g~')n. (16) 

In order to study the disorder-AF transition in the F4 lattice (/? < 0), 
we perform a change of variables which transforms the plane-AF vacuum 
(suppose X — y) into a ferromagnetic one defining 

Ky.t = i-lT'-'^.y.t. (17) 

Q matrix changes and then the propagator, Q~^, reads 
«'f> = Hy) ^- + 4(2 + ,(p)) - 



where 



g{p) = cospx cos Py + cospz cos pt — cosp^ cos Pz (19) 

— COSPx. cos Pt — COS Py COS Pz — COS Py COS Pt, 
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N 




1 


-0.17459(15) 


2 


-0.175608(5) 


3 


-0.17619(1) 


4 


-0.1766(1) 



Table 3: Critical couplings divided by for different O(A^) models. We 
include also the obtained in references [Q] and 



and ^ is defined from (translationally invariant) auxiliary fields as 

2(A, + i«,) = (-/5')(r' + 8). (20) 
From the variational equation ([T6|) we obtain f3'^ imposing ^ = oo: 



(3' = / —L- _ = -0.178972. (21) 

J (27r)4-8-4(7(p) ^ ^ 

In table |^ we compare ( pT]) with MC results for = 1, 2, 3, 4. The good 
agreement between the simulations for these values of and the analytical 
limit when iV — >■ oo points to the absence of an abrupt change of the critical 
properties as a function of A^. Exactly at A^ = cxd the order of the transition 
is not clear, because the divergence in the correlation length can simply be 
caused by the Goldstone bosons of the symmetry breaking. 



9 Conclusions 

In this letter we present a MC study of the four dimensional antiferromagnetic 
0(2) and 0(3) models in the F4 lattice. We study the critical behavior of 
these models with FSS techniques. There is an apparent asymptotic behavior 
which gives false critical exponents for not large enough lattice sizes. This 
transitory effect can be understood as caused by a large correlation length 
whose presence can be demonstrated for some observables (as the eigenvalues 
of the sum tensor of the period- two magnetizations). This must be very 
carefully controlled, because as we see in our case, the behavior changes 
drastically when larger lattice sizes are considered, revealing the true first 
order nature of the 0(2), 0(3) transitions. We also see that this effect is 
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bigger as N grows, so that for larger values of N, it is very difficult to 
study numerically the critical properties of the system. However, the great 
accuracy in the determination of the critical point obtained by the analytical 
calculation at N = oo points to a similar qualitative behavior for all values 
of N. 
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